Derealization of brane gravity by a bulk black hole 
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We investigate the analogue of the Randall-Sundrum brane-world in the case when the bulk 
contains a black hole. Instead of the static vacuum Minkowski brane of the RS model, we have 
an Einstein static vacuum brane. We find that the presence of the bulk black hole has a dramatic 
effect on the gravity that is felt by brane observers. In the RS model, the 5D graviton has a stable 
localized zero-mode that reproduces 4D gravity on the brane at low energies. With a bulk black 
hole, there is no such solution — gravity is delocalized by the 5D horizon. However, the brane 
does support a discrete spectrum of metastable massive bound states, or quasinormal modes, as was 
recently shown to be the case in the RS scenario. These states should dominate the high frequency 
component of the bulk gravity wave spectrum on a cosmological brane. We expect our results to 
generalize to any bulk spacetime containing a Killing horizon. 



Introduction 

The Randall-Sundrum (RS) model Q consists of a 4D 
Minkowski brane embedded (with mirror symmetry) in 
a 5D anti de Sitter (AdS) bulk spacetime. Although the 
extra dimension is infinite, it is exponentially warped and 
this leads to the recovery of 4D gravity on the brane at 
low energies. Thus the RS model provides the basis for 
models that could describe the observed universe as a 
brane-world. The key to this feature is the recovery of 
General Relativity at low energies, i.e., the existence of a 
normalizable zero-mass mode of the 5D graviton that is 
a bound state on the brane. General Relativity acquires 
corrections from the massive modes of the 5D graviton 
on the brane. (The massive modes dominate over the 
zero-mode at high energies.) 

Why does the RS model allow for a brane-localized 
zero-mode? The reason is the warp-factor, which erects 
a potential barrier around the brane that efficiently 
"squeezes" bulk gravity waves. In this work, we ex- 
plore the consequences of modifying this barrier by intro- 
ducing Weyl curvature into the model via a bulk black 
hole. In order to have a static vacuum configuration, 
we require a curved brane with an Einstein-static metric. 
The black hole's presence has a drastic implication for 
the zero-mode: there is no longer a gravity-wave bound 
state solution, because the extreme gravitational field 
near the horizon causes the brane potential barrier to 
become highly permeable, or "leaky" . 

The bulk metric satisfies the 5D Einstein equations 
Gab = ~ A 59a b \ and is given by 
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where Rq determines the ADM mass of the black hole, I 
is the AdS length scale defined by the bulk cosmological 



constant (Ag = — 6/£ 2 ), and 
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Thus, there is an event horizon at R = R^. We define 
dimensionless coordinates (t, r) = (T/R^, R/R^), so that 



m = 



(r 2 +7 2 + l)(r 2 - 1) 
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The bulk geometry is completely characterized by the 
ratio of the AdS length scale to the black hole horizon 
radius, 7, and the horizon is always at r = 1. It is sen- 
sible to call solutions with 7 < 1 "big" black holes and 
solutions with 7> 1 "small" black holes. 

A static brane is introduced by identifying r = as 
the boundary of the bulk and discarding the region r > 
tv The brane metric has the Einstein static form, 



(5) 



where r = fhRht = fhT is the cosmic time. In the RS 
model, the tension (vacuum energy) of the brane exactly 
cancels the bulk cosmological constant so that the ef- 
fective 4D cosmological constant vanishes and the brane 
is Minkowski. Here, the tension exceeds the critical RS 
value, so that the brane has positive cosmological con- 
stant. The acceleration due to this effective cosmological 
constant is nullified by the "dark radiation" induced by 
the projection of the bulk Weyl curvature onto the brane. 
Hence, unlike General Relativity we do not need any ad- 
ditional matter to have an Einstein-static configuration 
0- This is a delicate balance, which is why the brane's 
position must be fine-tuned. 

We can verify this via the junction conditions. The 
jump in the extrinsic curvature is determined by the en- 
ergy and stresses on the brane. By the standard junction 
conditions and the mirror symmetry, the extrinsic curva- 
ture of a brane with tension a must satisfy: 
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where n a = — f~ 1 ^ 2 6 a Ii is the normal. Together with 
Eq. (J5J, this leads to two equations in the three parame- 
ters 7, a, and r*b, with solutions 
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The solution for rb shows that a pure tension brane is 
coincident with the photon-sphere of the bulk black hole. 
(This is a special case of the general result that pure 
tension branes are always totally geodesic for null paths, 
i.e., they are umbilical surfaces.) Equations JJJ, (JSJ) and 
define the Einstein-static (ES) brane-world model. 



Tensor perturbations 



We now consider tensor perturbations of the bulk met- 



ric 
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the metric perturbation, 
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The tensor harmonics are defined by 
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where V,- is the covariant derivative of h 
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and 



k 2 = L{L + 2)-2, L = 1,2,3,... (11) 

Using the general results of Kodama and Ishibashi Q 
specialized to 5D Schwarzschild-AdS, we see that the lin- 
earized Einstein equations show that ip k satisfies a wave 
equation, 
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where x is the tortoise coordinate and the bulk black hole 
horizon is at x — — oo. The boundary condition follows 
from the requirement that the bulk fluctuation preserves 
the matter content of the brane, i.e., the brane has only 
tension after perturbation: SK a (, = —n 2 a 5g a b/6, which 
implies 



d r {r- z,2 ^ k ) 
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It is useful to compare Eq. (|12f> with the analogous 
master wave equation in the one-brane RS model: 

d 2 cj) k d 2 (f) k , 

U k {z)(t> k , (15) 



<9i 2 



dz 2 



fc 2 



15 



4(z-2) 2 



(16) 



Randall- Sundrum : 



15 



4(z-2) 2 



'1 



Einstein-static: 



Vfc oc exp(2«;a;) 



\ 



-3S(z - 1) 



-3/ b r b 1 ^(as- - ccb) 



FIG. 1: Comparison of the tensor potentials in the Randall- 
Sundrum and Einstein-static brane-worlds. 



where z is a dimensionless conformal coordinate. (For 
ease of comparison, we are using the mirror image of the 
bulk on the left of the brane, z b = 1.) Since the spatial 
geometry is flat in the RS scenario, k is a continuous 
parameter. The RS boundary condition is 



d z (z- 3/2 M 



= 0. 
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The U k and V k potentials are sketched in Fig. where 
we have added a delta function at the brane position to 
enforce the boundary condition. The major difference 
between the two potentials is their asymptotic behaviour 
far from the branes. In the RS case, the potential decays 
like 1/z 2 to a positive constant. By contrast, the ES 
potential decays as e 2KX for x — > — oo (the black hole 
horizon), where k is the dimensionless surface gravity: 



Derealization of the zero-mode 

The asymptotic behaviour of the potentials is crucial 
for the zero-mode. To see this, we assume e*"' time- 
dependence for ip k and <p k , which converts Eqs. (|12|) and 
(|15|l into Schrodinger-type eigenvalue problems with en- 
ergy parameter E = uj 2 . From elementary wave me- 
chanics, it is plausible for the RS potential to support a 
positive-energy normalizable bound state, because there 
are both an attractive delta function and an asymptoti- 
cally positive potential. Indeed, such a bound state does 
exist with lu = ±k; this is the zero-mode responsible 
for reproducing General Relativity on the brane. By 
the same token, the fact that the asymptotic ES po- 
tential vanishes strongly means that it cannot support 



a normalizable bound state with 



> 0. There is 



Randall- Sundrum type zero-mode in the Einstein- static 
brane-world. The vanishing of the potential reflects the 
fact that the horizon is perfectly transparent to gravity 
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FIG. 2: Numerical scattering experiment: a Gaussian gravity 
wave pulse collides with an Einstein-static brane. 



waves, so it is fair to say that the derealization of brane 
gravity is entirely due to the bulk black hole. 1 

Is there a tachyonic instability with uj 2 < in the ES 
brane-world? This cannot be determined from inspection 
of the potential in Fig. ^ O ne way to answer this is to 
conduct 'scattering experiments' where the wave equa- 
tion H12JI is solved numerically 0. As initial data, we 
choose a Gaussian pulse moving towards the brane at 
t = 0. Figure [3 shows a rather clean scattering event 
where the pulse strikes the brane, is reflected, and then 
propagates to infinity. At late times, the brane geometry 
reverts to its background configuration as the fluctuation 
dies away. This suggests that there is no linear instability 
in the system, and that none of the energy in the origi- 
nal pulse becomes trapped on the brane. Experimenting 
with a wide variety of incident signals suggests that both 
of these conclusions are generic and independent of initial 
data. 

More inferences can be drawn from the late-time grav- 
ity wave signals on the brane in these scattering exper- 
iments, which are shown in Fig. |3 We consider Gaus- 
sian incident pulses with identical characteristics in both 
the RS and ES (with k = L = 1) scenarios, and di- 
rectly compare the resulting waveforms. The RS signal 
is very well approximated by 4>k ~ Re(^4e lwt ), where 
A is a complex constant. A least-squares fitting gives 
lu = 1. 00, which implies that we are seeing stable zero- 
mode oscillations with u> = k. On the other hand, the 
ES signal is exponentially damped, in agreement with 
the behaviour seen in Fig. [2] The waveform is again 
well-described by ipk ~ Re(Ae luJt ), but with a complex 
frequency uj — 0.848 + 0.256 i. Hence as t — > oo, some of 



1 This conclusion is reinforced by the consideration of an ES brane 
embedded in pure anti-deSitter space, without a black hole. In 
that case, the 5D scalar fluctuation spectrum includes a stable 
brane-localized mode Q , in line with our claim that the horizon 
is responsible for the derealization of brane gravity. 
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FIG. 3: Brane gravity wave signals induced by the scatter- 
ing of identical Gaussian pulses in the Randall-Sundrum and 
Einstein-static (7 = 5) scenarios, with k = 1. The values of 
u> are obtained by fitting a Re( J 4e*" t ) template to each wave- 
form. 



the energy in the incident pulse remains localized on the 
RS brane while the ES brane radiates it all away. This is 
a direct numerical confirmation of the claim made above: 
the brane in the RS scenario supports a stable zero-mode, 
while the brane in the ES scenario does not. 



Quasinormal modes 

The removal of the zero-mode by the bulk black hole is 
our main result, but the complex- frequency oscillations 
exhibited in the ES brane-world deserve some investi- 
gation. Such behaviour is reminiscent of the familiar 
'ringdown' waveform from black hole perturbation the- 
ory, which is a direct consequence of the existence of 
so-called 'quasinormal modes' (QNMs). These are so- 
lutions of the relevant master wave equation subject to 
purely outgoing boundary conditions, and are charac- 
teristic of systems where energy can be lost to infinity. 
QNMs are described by a discrete set of complex frequen- 
cies u) n with Imw„ > 0. Hence, QNMs are exponentially 
damped in time. These modes are naturally interpreted 
as metastable bound states or scattering resonances of 
the potential, i.e., 'almost trapped' waves |6j. 

It has recently been demonstrated that the brane in 
the RS scenario supports QNMs Q, which are discrete 
modes embedded within the Kaluza-Klein continuum of 
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massive modes. Thus it is perhaps not surprising that 
the ES brane exhibits quasinormal ringing. Indeed, one 
might have expected QNMs from the fact that we are 
really considering a kind of modified 5D black hole per- 
turbation problem with a boundary. The novel feature 
is that the ES late-time gravity wave signals are dom- 
inated by QNM oscillations. This is in contrast to the 
RS scenario, where the zero-mode usually obscures the 
quasinormal ringing. Since these scattering resonances 
are of crucial importance to actual gravity wave signals, 
we calculate the QNM spectrum for the ES brane- world. 

Our method for finding the QNMs, which is discussed 
in detail elsewhere Q, relies on the series solution of 
Eq. I|12|l in the frequency domain. Assuming such a solu- 
tion satisfies both the brane boundary condition Eq. 1141) 
and the outgoing wave condition, 



4>k 



e Mt+x) ag x 
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results in an infinite-order polynomial in to. Truncation 
of this polynomial at some order TV gives a finite number 
of complex solutions uJn . Frequencies that are stable in 
the N — > oo limit are the QNM frequencies of the system. 

We apply this method to calculate the first eight quasi- 
normal frequencies of the ES brane in the L — 1 case, as 
functions of 7, and the results are plotted in Fig.0] Fre- 
quencies are labelled in order of increasing modulus, the 
smallest being the fundamental mode (n = 0) and the 
others being the overtones (n — 1, . . . , 7). We can check 
the validity of the results by examining the value of the 
fundamental frequency for 7 = 5: 



uj Q = 0.8488226912 + 0.2563765163 i, 



(20) 



which is in excellent agreement with the best-fit fre- 
quency in the bottom panel of Fig. ^ 

Figure 0] shows that the QNM frequencies approach 
constant values as 7 becomes large (i.e. for small black 
holes) . This can be understood by examining the leading 
order behaviour as 7 — > 00 of the master equation l|12|) 
and boundary condition i|14fl . In this limit, r < r\, w 
and we can neglect terms of order r 2 /j 2 . This results in 
the asymptotic wave equation, 



d 2 j> k = 
dt 2 

where g(r) 
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1 — 1/r 2 . The boundary condition reduces 



dr 



2V2 



= 0. 



(22) 



Both equations are independent of 7, and thus the QNM 
frequencies should also be independent of 7 in the 7 3> 
1 limit, as confirmed in Fig. We also find that the 
frequencies are evenly spaced for large 7, and are well 
approximated by: 



uj n w 0.78 + 0.28 i+ (0.57 +1.81 i)n. 



(23) 
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FIG. 4: Real (top) and imaginary (bottom) parts of the first 
eight quasinormal frequencies in the ES brane-world. 



The root mean square error between this relationship and 
the calculated frequencies is 0.02 at 7 = 100. 

In the small 7 limit, since r > 1, we can neglect terms 



of order 7 2 /r 2 . Writing ipk(t,r) 
equation for 7 <C 1 is 



e **(r), the wave 



-7 4 w 2 * fc = 



dr 



/i(15r 4 



9) 



4r 4 
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where h(r) = (r — l)/r 2 . In this equation, 7 2 and w 
are degenerate because they only appear in the product 
-f 2 oj. This suggests that the QNM frequencies of the sys- 
tem should obey a power-law scaling u> n oc 7 -2 ss k for 
7 <C 1. This is indeed true for the case of scalar Q, elec- 
tromagnetic, and gravitational || fluctuations of a pure 
Schwarzschild-AdS black hole with no brane around it. 
However, the brane boundary condition for 7 <C 1, 



37 
2V2 



*fe(r b ) 



breaks the degeneracy, and we should not necessarily ex- 
pect that uj n cx 7~ 2 for 7^0. But in fact, we do see 
such a scaling for the overtone frequencies in Fig. 



n nl - 2 for 7 < 0.14, n = 1,2, 



(26) 



where VL n is a complex constant determined numerically. 
For example, by performing a fit between 7 = 0.10 and 
0.14, we find Cli = 3.109 + 2.814 i. The goodness of the fit 
can be assessed by looking at the RMS discrepancy be- 
tween the logarithms of the calculated and approximate 
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Kaluza-Klein masses? 



FIG. 5: The first ten quasinormal frequencies of the 7 = 5 ES 
brane for various values of L. 



In the RS scenario, the various modes of the 5D gravi- 
ton can be labelled by an effective 4D Kaluza-Klein (KK) 
mass. This can be understood as follows: The RS poten- 
tial in HI 5(1 can be re-written as 



U k {z) = k 2 + U{z) 



(28) 



where U is independent of the spatial scale k. With 
4>k(t, z) = e iult ipk(z), we find the dispersion relation 



2 ci2 2 

m = a — p , 



E = lj/£, p = k/£, 



(29) 



where m 2 £ 2 is the eigenvalue of —(d/dz) 2 + U(z). To a 
stationary brane observer, E is the energy of the mode 
while p is its 3- moment urn. Hence, m = y 1 E 2 — p 2 is 
naturally interpreted as the effective 4D graviton mass 
according to the standard KK paradigm. 

By contrast, in the ES brane-world, the spatial scale k 
can no longer be separated from the potential Vk, and no 
simple position-independent dispersion relation like 1(29(1 
can be found. At best, we can define a local dispersion 
relation on the brane: 



frequencies. For n — 1, this error is 10 4 and 2 x 10 3 
for the real and imaginary parts, respectively. 

The small-7 behaviour of the fundamental mode is 
quite different. We find 



~ 2 7Tl2 ~2 

m = E — p , 



E = 



-RhVJb 



P = 



Rhrb 



(30) 



ujq w I.687- 1 + 0.497 % for 7 < 0.14, 
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i.e., the real part appears to scale like 7 _1 while the imag- 
inary part approaches a constant. The discrepancy in 
the asymptotic behaviour of the overtones i(26|) and the 
fundamental mode l(27|) would seem to suggest that the 
latter is more sensitive to the boundary condition. 

Attempts to test these relations for much smaller val- 
ues of 7 are constrained by computing speed, since for 
7 — * 0, we have n> — ► 00, which is a singular point of 
the master wave equation 1 (12(1 in the frequency domain. 
Hence the series solution for ipk becomes poorly conver- 
gent at the brane, which means we must retain unrea- 
sonably many terms to get accurate QNM frequencies. 

We have also calculated QNM frequencies for different 
values of L, i.e., for perturbations on different spatial 
scales on the brane. The results for 7 = 5 are shown 
in Fig. [5] The general trend is that the real parts of the 
frequencies increase with L (i.e., for smaller scales), while 
the imaginary parts remain roughly constant or decrease 
slightly. This is in keeping with results obtained for the 
scalar QNMs of a 4D Schwarzschild-AdS black hole with 
no brane [8j . 



Note that in the shortwave approximation k 3> 1, E and 
p are the mode's energy and 3-momentum as measured 
by comoving brane observers, respectively, and m is the 
effective mass of the mode. 

What are the masses of the QNMs we calculated 
above? We continue (j3*0"|) into the complex plane in order 
to define a complex mass for each QNM. Then, Re?7i rl 
is interpreted as the conventional mass of the resonance, 
while (Imm„) _1 is roughly its half-life. Our results for 
m n for L = 80 are shown in Fig. El Interestingly, the up- 
per panel suggests a mass gap between the fundamental 
mode and the overtones for small 7, i.e., for large black 
hole effects. This is reminiscent of a de Sitter brane- 
world [13 (although in the de Sitter case, the gap is be- 
tween the zero-mode and the KK continuum). 



Generalizations 

One can generalize the ES brane-world model by mak- 
ing the bulk more complicated, letting the brane move, or 
both. Will our results carry over to these new situations? 

First consider the addition of matter fields in the bulk. 
These may take the form of dilatons, moduli fields, su- 
pergravity form fields, etc. Stationary, spherically sym- 
metric solutions sourced by such matter and featuring 
an event horizon are classified as "dirty black holes" . 
Scalar wave propagation in 4D dirty black hole space- 
times shows j 1 11 that the potential in the master wave 
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1.6i n = 7 
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FIG. 6: Effective 4D masses, as measured by brane observers, 
of the first eight QNMs in the L = 80 case 

equation vanishes at the horizon, just as in Eq. (|12J) . 
Now, the key reason that the zero-mode becomes delocal- 
ized in the ES model is the fact that V(Rh) = 0. Hence, 
if we make the natural assumption that the behaviour 
seen in Ref. generalizes to spin-2 fields and higher- 
dimensional backgrounds, we conclude that the zero- 
mode will remain delocalized if the 5D Schwarzschild- 
AdS bulk is replaced with a dirty black hole. Indeed, 
we can go even further by conjecturing that any static 
brane located outside a stationary Killing horizon can- 
not support a normalizable zero-mode, precisely because 
any such horizon be will completely transparent to bulk 
gravity waves. 

Next, consider brane motion in a Schwarzschild-AdS 
bulk, rb = r\){t) as in model of brane cosmology. Do 
the QNM frequencies we have calculated tell us any- 
thing about the behaviour of cosmological perturbations? 
Fluctuations with |w| ^> H = fb/^b do not 'feel' the ex- 
pansion of the universe and effectively 'see' the brane as 
stationary, i.e., the characteristic timescale of the pertur- 
bation is much larger than the characteristic timescale of 
the cosmological dynamics. Hence, we can expect the 
QNMs calculated for a static brane to be approximate 
solutions for the gravity waves around a moving brane 
if l/|w„| is much shorter than the Hubble time. Such 
QNMs will dominate that part of the bulk gravity wave 
spectrum. Hence, the high frequency lo H component 
of the bulk gravity wave spectrum on a cosmological brane 
should be dominated by the metastable bound state reso- 
nances of the corresponding static brane. 

Such an approximation has the potential to greatly 
simplify the thorny problem of brane cosmological per- 



turbations, but there is a significant caveat. The QNMs 
we have calculated are only for one brane position, i.e., on 
the photon sphere. In order to have a complete picture, 
we need to know the QNM frequencies for static branes 
over a range in r, corresponding to the addition of matter 
on the ES brane. The calculation of these frequencies is 
the subject of a separate paper 0. 



Conclusions 

We have considered static pure-tension branes sur- 
rounding a 5D bulk black hole. By studying the tensor 
perturbations, we have seen that the bulk Killing hori- 
zon causes the brane's zero-mode to become delocalized. 
In other words, the gravitational field of the black hole 
makes it impossible for the brane to support a normaliz- 
able bound state. However, we also found that the brane 
supports a discrete spectrum of metastable bound states, 
or quasinormal modes, as in the Randall-Sundrum sce- 
nario. Using a series solution of the master wave equa- 
tion, we have calculated the quasinormal frequency spec- 
trum. We discussed why the massive mode Kaluza-Klein 
decomposition common in other brane- world models does 
not work in the current problem, but then showed how 
one could define an effective local mass measured by 
brane observers. The locally defined mass shows a gap 
between the fundamental and overtone modes. 

Our results are expected to generalize in several impor- 
tant ways. We expect that whenever there is a stationary 
Killing horizon in the bulk, a surrounding brane cannot 
support a normalizable bound state. Furthermore, we 
expect that the high frequency bulk gravity wave spec- 
trum on a moving brane will be well represented by a 
sum over the quasinormal resonances of the correspond- 
ing static brane. 
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